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IM'TRODUCTION 


1 1 


The  purpose  of  this  r ludy  is  to  determine  the  frequency  spread 
caused  by  reflection  of  an  acoustic  signal  off  a time-varying  stochastic 
.'■urface.^  Scattering  from  a rough  surface  has  been  investigated  by  many 
scientists,  of  which  we  will  mention  only  a few.  Eckart,  using  the 
method  of  Physical  Optics  and  applying  the  Helmholtz  integral,  obtained 
the  intensity  of  the  scattered  sound;  from  this  quantity,  the  scattering 
coefficient  of  the  surface  was  obtained.  Although  Eckart  had  the  time 
variation  of  the  surface  in  his  preliminary  integral,  he  did  not  com.pute 
any  quantities  involving  the  time  variation.  Beckmann,  in  his  book, 
also  applied  the  Heliriholtz  integral  to  a rough  stochastic  surface,  but 
did  not  consider  a tim.e-varying  surface.  Marsh,  on  the  other  hand, 
utilized  the  Rayleigh  mef.hod  and  considered  the  scattered  wave  to  consist 
of  an  infinite  set  of  outgoing  plane  waves.  By  applying  boundary 
conditidns  at  the  surface,  the  scattering  coefficient  of  the  surface 
was  obtained. 

The  number  of  investigations  of  a time-varying  surface  is  consider- 
ably less  than  that  of  a fixed  surface.  Of  these,  we  will  mention  only 
the  works  of  Roderick  and  Parkins.  Roderick  considered  a sinusoidal 
traveling  surface  wave.  For  a single-frequency  incident  acoustic  wave, 
the  frequency  of  the  reflection  in  the  specular  direction  i“"  the  same 
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cs  the  incident  frequency,  end  the  frcquencier  in  the  non-npeculnr 
directions  nre  equal  to  the  incident  frequency  plus  a multiple  of  the 
surface  frequency.  Parkins  considered  a random  time-varying  .surface, 
a joint  Gaussian  dis tritu tion  of  heights,  and  a .’Jeumann-l  ierson  surface 
height  spectrum.  For  both  small  and  large  surface  roughness,  he  obtained 
the  scattered  acoustic  spectnam. 


In  our  work,  we  start  with  the  assumptions  used  by  others:  the 

method  of  Physical  Optics  is  used;  the  insonified  area  is  in  the  far 
field  of  the  source;  the  observation  point  is  in  the  far  field  of  the 
inso.nified  area;  the  surl'ace  is  pressure  release;  no  shadowing  or 
multiple  reflections  occur:  a directional  source  is  used;  and  variation 
of  the  surface  is  slow  compared  to  the  speed  of  sound. 


(Figure  l)  Starting  with  the  Helmholtz  integral,  and  a single- 
frequency acoustic  source,  the  com.plex  press-ire  at  the  observation  point 
is  obtained.  It  is  assumed  that  the  surface  is  randomi,  homogeneous,  and 
stationary.  That  is,  for  a given  instant  of  time,  first-  and  second-order 
properties  of  the  s-urface  are  independent  of  absolute  position,  and  depend 
only  on  the  differences  of  the  x and  y coordinates,  these  differences 
being  specified  by  u and  v.  In  time,  the  properties  depend  only  oti  the 
difference  in  tim.e,  specified  by  T • The  autocorrelation  of  the 
received  complex  pressure  is  obtained  and  is  denoted  by  Rp(T^.  C depends 
on  the  geomietry  of  the  situation,  a,  b and  c are  the  sums  of  the  direction 
cosines  in  the  x,  y and  z directions,  respectively,  of  the  incident  and 
scattered  waves. 

is  the  acoustic  wave  number. 


0~  is  the  rms  wave  height. 


(u,v,  t)  is  the  normalized  surface  correlation  f'unction, 
be  discussed  on  the  next  slide. 


which 


P(t)  is  the  expected  value  of  the  received  com.plex  pressure, 
and  is  re; erred  to  as  the  coherent  component  of  the  received  acoustic 
pressure  waveformi.  We  have  assumed  a joint  Gaussian  distribution  of 
heights  in  this  equation. 


a 


The  Fourier  transform  of  Rp(Tn  gives  the  received  acoustic  power 
spectrum.  The  Fourier  transform  of  the  first  term  ol’  results  in 

a line  component  located  at  the  transm.itted  frequency.  This  line 
component  has  appreciable  value  mainly  in  the  specular  direction  and 
only  for  a slightly  rough  surface.  For  a given  geometry  and  given 
surface  correlation  function,  the  doul'le  integral  in  the  slide  could  be 
evaluated  by  numerical  methods. 
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(p'i(*ure  2)  ^{^,^1'^')  ropr'^ront.r;  the  s'lrf’iee  wive  hoLi'tit  -it  the 

point,  at  time  -t  , anfi  |C7(u,v,  T)  in  the  nortnrn  Lzeh  npatia] -temporal 

correlation  of  the  surface  wave  heif'hts.  As  stated  ear.ier,  ! arkins 
has  obtained  the  receiveil  a.coustic  power  spectrum  for  the  care  of  the 
h'eumann-Pierson  surface  height  spectrum  and  a fui  ly-aroused  sea.  This 
corresponds  to  a wide-hand  spectrum,  with  a approximately  ecuai  to 
1.4.  We,  instead,  consider  a narrow-band  surface  spectr:im,  which 
corresponds  to  the  case  of  swell  waves.  For  the  narrow-band  ass’imption, 
the  correlation  of  surface  heit'hts  is  in  the  form  {^iven  by  the  second 
equation.  The  magnitude  and  phase  of  t.he  complex  envelope  /9(u..v,t) 
vary  slowly  in  ~C  with  respect  to  cos  ( 7-rrf^T  ).  If,  in  addition, 
the  iso-values  of  the  correlation  function  is  an  ell  ipse,  the  phase  of 
^ (u,v,0)  is  identically  zero.  In  the  computations  to  be  performed, 

;n  ellipticcal  iso-correlation  surface  is  assumed. 


(Fig’ure  3)  In  the  expression  for  the  received  power  spectrixm, 
there  is  a term  that  contains  an  exponential  of  the  correlation  function 
sVio'wn  in  the  previous  slide.  For  the  narrow-band  case,  this  correlation 
is  approximately  sinusoidal,  and  the  expression  can  t>e  expanded  in  a 
series  of  Bessel  functions  of  imaginary  argument.  An  examination  of 
the  resulting  expression  shows  that  the  received  acoustic  spectrum 
consists  of  a series  of  spectral  lobes  or  sidebands,  each  separated 
from  the  incident  acoustic  frequency  by  multiples  of  the  center  frequency 
of  the  surface  variation  plus  a spike  at  the  acoustic  frequency  ^ . 
The  coherent  component  occurs  mainly  in  the  specul.ar  direction  and  only 
has  an  appreciable  value  for  small  surface  roughness.  We  will  now  find 
the  power  in  each  sideband  component. 


(Figure  4)  The  relative  scattered  power  in  each  sideband  can  be 
obtained  in  terms  of  only  two  parameters  of  the  surface,  namely,  and 

is  the  usual  measure  of  the  roughness  of  the  surface  . 

Ly  and  L are  the  (correlation)  distances  along  the  x and  y axes, 
at  which  the  correlation  of  surface  heights  is  i/e  of  its  ps.ak  value. 

^ in  this  equation  represents  the  spatial  dependence  of  the  complex 
envelope  of  the  surface  correlation  function. 


S is  the  Kronecker  delta. 


gives  the  relative  scattered  power  in  t 
m = 0 corresponds  to  the  lobe  centered  at  acoustic 
scattering  coefficient  is  defined  by  Kckart, 
factor  C*,  as  given  by  Hort.on  and  Muir,  are  shown 
equation.  This  equation  contains  the  factors  from 


he  m^^^  sideband, 
frequency  f,^.  The 
and  the  correction 
iti  the  second 
which  the  scattering 
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r.trerif-Tth  couH  be  evaluated  for  a rpeoific  (geometry.  B in  a factor 
which  dependn  only  on  geometry. 

t^(0)  is  not  integrable  in  closed  form.  A computer  program  has 
been  ’writ, ten  to  ntunerically  evaluate  ^(O)  for  given  values  of  ^ and 
^ . An  investigation  into  the  accuracy  of  the  numerical  integration 

iy  Simpson's  rule  shows  that  from  6 to  8 places  of  accuracy  are  obtained. 
The  number  of  summations  for  each  value  of  o<  and  ^ ranged  from  400  to 
4,000.  '/ith  the  advent  of  modern  com.puters,  this  tpomerical  integration 
is  feasible  and  economical. 


(Figure  5)  In  our  computations,  we  consider  two  types  of  auto- 
correlations of  surface  height,  n-'imely,  exponential-  and  Gaussian- 
modulation  of  an  oscillatory  function.  In  this  fi^^e,  we  have  exar.ined 
G" ussian-raodulation  of  a sinusoid  for  o<  equal  to  zero;  this  corresponds 
to  a and  ! , the  direction  cosines  with  respect  to  the  x and  y axes, 
equal  to  zero,  and  thus  to  the  specular  direction.  The  power  in  the 
first  sideband  is  plotted  as  a function  of  p . ^-O  corresponds  to  a 
flat  surface;  for  this  value  of  |9  , there  is  no  scattered  power  in  the 
sidebands,  as  all  the  power  appears  in  the  coherent  com.ponent.  As  A 
increases,  ther'c  is  less  power  in  the  coherent  component  and  more 
scattered  power  in  the  first  sideband.  The  power  in  the  first  sideband 
reaches  a maxim.um  value  when  the  surface  roughness  parameter  S is 
approximately  unity,  and  tiien  decreases.  For  these  larger  vaWes  of  p , 
the  surface  is  rougher,  and  more  of  the  scattered  power  will  appear  in 
higher  order  sidei>ands. 


(Figurf’  6)  In  this  final  figure,  a three-dimensional  plot  of  the 

i scattered  pov;*' r in  the  first  sideband  is  drawn  vs  and  p'  . For  any 

value  of  eg  , as  increases,  the  power  in  a particular  sideband  will 
• first  increase  and  then  decrease.  For  a given  value  of  p'  , i.e. 

i roughness,  as  ot  increases,  again  the  power  first  increases  and  then 

decreases.  Three-dimensional  plots  of  this  type  were  obtained  for  several 
sidebands  and  for  Loth  Gaussian-  and  exponentially-modulated  oscillatory 
functions. 


Before  the  final  conclusion  and  summary,  we  would  like  to  discuss 
the  difference  between  the  narrow-band  Gaussian  surface  and  the  fixed- 
amplitude  sinusoidal  surface  considered  by  Roderick.  We  have  considered 
a narrow-band  spectrum  for  the  height  variation  at  a point  of  the  surface, 
and  assumed  ♦he  joint  probability  dehsity  of  the  surface  heights  to  be 
Gaussian.  As  the  bandwidth  of  the  surface  variation  decreases  and 
approaches  zero,  the  properties  of  this  surface  process  do  not  approach 
tlie  fixed-amplitude  sinusoidal  surface  case.  For  example,  the 
distribution  of  heights  remains  Gaussian  and  is,  therefore,  different 
from  the  probability  density  associated  with  a sinusoid.  The  scattered 
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powern  in  the  rL3f!'in3r-  en  piven  by  the  narrow-band  Gaunsian  theory  are 
dll'Torent  from  thone  of  the  s iri{^le-frequency  ninu'oidal  theory. 

In  mathem-i ♦ ' vil  terms,  the  first-order  characteristic  fimction  of 
a f ixed-ampl it’jie  rinnsoid  is  a Bessel  function.  If  we  averat'e  this 
Bessel  function  with  the  Rayleigh  distribution  of  sine  wave  amplitudes, 
w^‘  oitair;  the  biussian  first-order  characteristic  function.  Thus,  our 
results  can  be  viewed  as  an  average  over  the  sinusoidal  surface  case. 

On  the  otlier  hand,  a narrow-band  time  function  with  center 
frequency  f^  and  : aniwidth  'A  resembles  a sine  wave  of  frequency  fg 
over  a short  perLs-d  of  time,  but  will  differ  from  a sine  wave  over  a 
lonper  period  of  time,  the  difference  showing  up  on  the  order  of  W sec. 
Thun,  if  we  ol  t-,in  a -h  rt  sample  of  the  reflected  signal  from  a narrow- 
band  surface,  an  i lotermine  the  received  acoustic  spectr'o.m  from  it,  the 
single-frequency  theory  (with  the  current  surface  height  amplitude) 
should  accurately  predict  that  spectram.  However,  for  a long  time 
average  over  the  statistical  properties  of  the  surface,  the  narrow-band 
Gaussian  theory,  not  ♦he  fixed -amplitude  single-frequency  theory,  must 
be  used  to  predic*  the  expected  behavior  of  the  spectruni. 

In  surLTiary,  we  have  considered  a time-varying  stochastic  surface, 
stationary  in  both  space  and  time.  The  reflected  acoustic  spectr’um 
was  obtained  in  terms  of  the  first-  and  second-order  characteristic 
f'unctions  of  the  surface  height  variations.  For  the  special  case  of 
surface  wave  height  variation  possessing  a narrow-band  spectrum  at  a 
point,  joint  Gaussian  distributions  of  heights,  an  elliptical  spatial 
correlation,  the  integral  expression  obtained  for  the  reflected  acoustic 
sidelobe  powers  at  the  observation  point  was  derived  and  evaluated.  The 
integral  was  computed  in  terms  of  two  fundamental  parameters,  and  curves 
of  the  results  were  shown. 
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SCATTERED  POWER  IN  THE 


SCATTERED  POWER  IN  THE  SPECULAR  DIRECTION 
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SCATTERED  POWER  IN  THE  FIRST  SIDEBAND 


o 

fc-: 

0) 

M 

) 

C-  r— 4 

<D 

E-  <r\ 

I 

r-J 
: ■}  r-A 
• ^ Tsi 
■ ■: 


Navy  Underwater  Sound  Laboratory  Official  Photograph 

• 37350  - 10  - 69 


